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ABSTRACT 

Numerical simulations suggest that active galactic nuclei (AGNs) play an important role in the formation 
of early-type galaxies by expelling gas and dust in powerful galactic winds and quenching star formation. 
However, the existence of AGN feedback capable of halting galaxy-wide star formation has yet to be obser- 
vationally confirmed. To investigate this question, we have obtained spectra of 14 post-starburst galaxies at 
z~0.6 to search for evidence of galactic winds. In 10/14 galaxies we detect Mg II AA2796,2803 absorption 
lines which are blueshifted by 490 - 2020 km s~'with respect to the stars. The median blueshift is 1 140 km s" 1 . 
We hypothesize that the outflowing gas represents a fossil galactic wind launched near the peak of the galaxy's 
activity, a few 100 Myr ago. The velocities we measure are intermediate between those of luminous starbursts 
and broad absorption line quasars, which suggests that feedback from an AGN may have played a role in 
expelling cool gas and shutting down star formation. 

Subject headings: galaxies: evolution — galaxies: ISM — galaxies: starburst — quasars: absorption lines 
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1. INTRODUCTION 

There is mounting evidence linking quasar activity to 
merger-induced star formation, but the precise timing and the 
physical relationship between the two are not well understood 
(Canalizo et al. 2006). Numerical simulations suggest that 
mergers of gas-rich galaxies induce radial gas inflows which 
fuel central star formation and black hole accretion. Sub- 
sequently, feedback from the active galactic nucleus (AGN) 
removes the gas and dus t and quenches star f ormation and 
black hole activity (e.g., iDi Matteo et al.ll2005l) . Such mod- 
els have enjoyed great popularity due to their success in re- 
producing the present-day properties of early-type galaxies, 
notably the color-magnitude relation and the correlation be- 
twee n black hole mass and bulge stellar veloci t y dispersion 
(e^ iGranato et all 120041: ISpringel etail 120051: Ovlenci et all 
2006). However, the existence of AGN feedback capable of 
halting galaxy-wide star formation has yet to be observation- 
ally confirmed. 

AGN feedback is predicted to quench star formation by 
re-heating the cold gas and expelling much of it in power- 
ful galactic winds. Galactic winds with velocities of 50 - 
500 km s" 1 are commonly detected in starburst galaxies via 
the presence of gas absorptio n lines that are blues hifted rel- 
ative to stellar features (e.g.. iHeckman et afluOOOl) . AGN- 
driven winds are expected to produce similar observationa l 
signatures, but higher outflow velocities dThacker et al.l2006h . 
The maximum feedback impulse is predicted to occur during 
the bright quasar phase. However, at this stage the quasar 
outshines the host galaxy, and provides a more ambiguous 
probe of the galaxy's interstellar medium (ISM). (A parsec- 
scale cloud near the quasar would be indistinguishable from 
a kiloparsec-scale galactic wind.) We have therefore elected 
to look for remnants of AGN-driven galactic winds during the 
post-starburst phase, a few 100 Myr after the peak of the star 
formation and AGN activity. 

1 Observations reported here were obtained at the MMT Observatory, a 
joint facility of the University of Arizona and the Smithsonian Institution. 

2 Steward Observatory, 933 N. Cherry Ave., Tucson, AZ 85721; 
tremonti @ as . arizona.edu, adiamond @ as . arizona.edu 

3 Hubble Fellow 

4 Center for Cosmology and Particle Physics, New York University, New 
York, NY 10003; jmoustakas@cosmo.nyu.edu 



Post-starburst galaxies are characterized by strong stellar 
Balmer absorption from A-stars, but little nebular emission 
indicative of on-going star formation. Local post-starbursts 
(sometimes called 'E+A' or 'K+A' galaxies) have the kine- 
matic and morphological signatures of pressure-supported 
spheroids, but frequently exhibit low surface brightness tidal 
tails i ndicative of a rece n t major merger (Zabludo ff et al.1 
Il996t iNortonet all 120011: I Yang et all |2004|) a nd signs of 
weak AGN activity (lYan et al.ll2006t lYang et al. 2006). Post- 
starbursts are therefore presumed to be late-stage mergers that 
have passed through their quasar phase and are in transition 
to becoming early-type galaxies. As such they provide ideal 
testing grounds for AGN-feedback models. 

We have obtained spectra of 14 post-starbursts at z ~ 0.6 in 
order to search for evidence of galactic winds that may have 
played a role in shutting down star formation. We have se- 
lected galaxies at intermediate redshift because this may be an 
important epoch f or the formation of early-type galaxies (e.g., 
iFaber et al.ll2005l) . and because the rest-frame near-UV is ac- 
cessible in the optical. Coverage of the near-UV improves 
our ability to estimate the recent star formation history of our 
galaxies and it enables us to measure Mg II AA2796,2803 
which is a sensitive probe of the ISM. We describe our ob- 
servations in §2, the stellar populations of our post-starbursts 
in §3, and the outflow kinematics in §4. We discuss evidence 
that the wind is powered by an AGN in §5 and conclude in §6. 

2. OBSERVATIONS AND DATA REDUCTION 

Our sample was selected from the Sloan Digital Sky Sur - 
vey (SDSS) Data Release 4 dAdelman-McCarthv et a02006h . 
The parent sample is composed of i < 20.5 mag objects that 
were targeted for SDSS spectroscopy as quasar candidates, 
but which were subsequently classified as galaxies at z= 0.5 - 
1. The typical signal-to-noise (S/N) ratio of the SDSS spectra 
is rather poor (S/N^2 pixel" 1 ), but sufficient for us to select 
a sample of galaxies for follow-up. We selected 159 objects 
with post-starburst characteristics — strong stellar Balmer ab- 
sorption and weak nebular emission (Tremonti et al., in prep.). 

We obtained higher S/N spectra of 14 galaxies with the 
Blue Channel Spectrograph on the 6.5-m MMT. We used 
the 500 line mm" 1 grating blazed at 5600 A which gave us 
spectral coverage from 4050 - 7200 A with a dispersion of 
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FIG. 1 . — Example spectra (black) and continuum model fits (gray). In the lefthand panel, the continuum models are offset for clarity. The righthand panel 
highlights the region around the Mg II doublet. Dotted lines mark the rest wavelength of Mg II . The presence of blueshifted lines indicates an outflow. The 
velocity of the most blueshifted component is given in km s~'in the upper right corner. 



1.19 A per pixel. For our z = 0.51-0.75 galaxies, this yielded 
rest-frame coverage from 2700 - 4100 A. Because most of 
our objects are unresolved in the SDSS imaging we used 
a 1" slit, which yielded a FWHM resolution of ^3.6 A. 
The spectra were reduced, extracted, and spectrophotomet- 
rically calibrated using the ISPEC 2D data reduction pack- 
age dMoustakas & Kennicuttll2006l) . The MMT spectra show 
good agreement with the SDSS data, but have S/N = 15-30 
pixel" 1 . Spectra of four representative galaxies are shown in 
Figure 1. 

3. STELLAR POPULATION 

We detect significant Mg II AA2796,2803 absorption in all 
of our galaxies. Mg II is one of the strongest interstellar res- 
onance absorption lines; however, it is also present in the at- 
mospheres of A-stars and later spectral types. Hence to ac- 
curately measure the ISM absorption lines , we must carefully 
model the stellar continuum. We use the iBruzual & Chariot 
d2003l hereafter BC03) stellar population synthesis models to 
create synthetic spectra for different star formation histories. 
We assume super-solar metallicity (Z=2.5 Z©) since Oxygen 
abundances several times solar are measure d in comparably 
lumin ous star forming galaxies at z ~ 0.6 (iLamareille et al.1 
2006). We adopt a star formation history designed to emulate 
a major merger between gas-rich disk galaxies. Following the 
starburst, star formation decays exponentially with time con- 
stants ranging from Tk„ rs ,=25 - 500 Myr. We fit each of our 
spectra with a grid of models spanning a range of ages, Tb U rst> 
and reddening values and adopt the model with the minimum 
X 2 as the best fit. In several cases to achieve an optimal fit 
it was necessary to add an additional power-law component 
which may represent a featureless quasar continuum (see §5). 

A byproduct of our stellar continuum modeling is an es- 
timate of the stellar mass, the time since the peak star for- 
mation event (tburst), and how quickly star formation ceased 
(Tburst)- These parameters and their uncertainties will be dis- 



cussed fully in Tremonti et al., in prep. Our modeling suggests 
that the galaxies are massive (0.7-4.8 x 10 11 M©) and have 
recently experienced a burst (tburst = 75 - 300 Myr) that faded 
rapidly (Tburst = 25- 100 Myr). The short starburst timescales 
imply strong feedback from supernovae or an AGN. 

4. GAS KINEMATICS 

We use our best-fit synthetic spectra to correct for the 
contribution of stellar absorption to the Mg II lines. At 
wavelengths l ess than 3300 A, the BC03 models use the 
iPicklesi d!998l) stellar library which has a spectral resolution 
of 10 A. This resolution is too low to adequately model the 
Mg II doublet in our data. We circumvent this problem by 
patching the BC03 models in the 2600-3300 A range using 
theoretical stellar spectra from the UVBLUE stellar library 
dRodrfguez-Merino et al.|[2005b . The rightmost panel in Fig- 
ure 1 shows our continuum fits in the 2750-2850 A region. 
In many cases the ISM lines are so strong or blueshifted that 
the stellar component of Mg II is unimportant. However, in a 
few galaxies stellar Mg II is dominant. After correcting for the 
stellar light we find that 10 of our 14 galaxies have measurable 
Mg II absorption. The equivalent widths (EWs) of interstellar 
Mg II range from 0.8 - 10.4 A (see Table 1). 

After correcting for the stellar contri bution to Mg II , w e 
fit the ISM absorption lines following Rup ke etail d2005"5) . 
In the optically thin case the doublet ratio is 2:1, but in our 
data the lines are moderately saturated. At our spectral res- 
olution (~ 100 km s" 1 ) this produces degeneracies between 
the optical depth at line center, the covering factor, and the 
Doppler b parameter. However, velocities can be measured 
robustly. The lines have a median Doppler width of b = 
260 km s" 1 , although it is possible that they include narrower 
unresolved components. In six of the galaxies we fit two ab- 
sorption components. Three galaxies display a P-Cygni pro- 
file — blueshifted absorption coupled with redshifted emis- 
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sion which may originate on the back side of an expanding 
shell. We model the emission with a Gaussian. In Table 1 we 
list the measured absorption-line velocities. We denote the 
velocity of the most blueshifted component in each spectrum 
as v max - We compute the average velocity, v avg , weighting the 
components by their EWs. The median values for the sample 
are v avg = -920 km s"'and v„ Mt = -1 140 km s" 1 . 

We hypothesize that the blueshifted Mg II lines originate 
in galactic winds that were launched near the peak of the 
galaxies' starburst activity a few 100 Myr ago. An alter- 
nate interpretation is that the absorbing gas is tidal debris as- 
sociated with the merger. However, the inner parts of tidal 
tails are e xpected to be bound and to fall back within a few 
100 Myr (Hibbard & Mihosl [19951) . whereas we detect out- 
flows. In addition, gaseous tidal tails are confined to relatively 
thin streams with small global covering factors; therefore it 
seems unlikely that we would detect tidal gas in absorption in 
70% of our sources. Hence we conclude that the blueshifted 
Mg II lines originate in fossil galactic winds. 

5. DISCUSSION 

Our z ~ 0.6 post-starburst galaxies rank among the most 
luminous and massive galaxies in the universe (M^ = -22.5 
- -23.7 mag, M* = 0.7-4.8 x 10 11 M ; Table 1) and they 
offer a rare window on the formation of today's massive 
early-type galaxies. The detection of interstellar Mg II in 
10 of our 14 galaxies enables us to probe the properties 
of the cold ISM. We find evidence for strong outflows in 
all 10 systems, with velocities in the range \ max = 500 - 
2000 km s -1 . The median v max of the sample is 1 140 km s -1 , 
which exceeds the 400 - 600 km s" 1 velocities typical of lumi- 
nous starburst galaxies (Heck man et al.ll2~000l) . We put these 
outflows in context in Figure 2 where we plot absolute B- 
band magnitude versus ISM outflow ve locity for a variety 
of sy s tems. We include loc al starbursts dSchwartz & M artin 
l2004t [Schwartz et al.ll2006l) . Lumi nous and Ultra-Lum inous 
Infrared Galaxies (LIRGs/ULIRGs: IRupke et al.ll2005bl), an d 
z — 3 Lyman Break Galaxies (LBGs; iPettini et all 1200 ll) . 
We augment this sample with s tarburst/AGN composite 
ULIRGs from Rup ke et ail d2005d) and a sample of Low- 
ionization Broad Absorption Line quasars (LoBALs) from 
the SDSS dTrump et alJ l2006). LoBALs are characterized 
by broad Mg II absorption troughs. They are more com- 
mon in infrared-selected than optically-selected quasar sam- 
ples dBoroson & MeyersH l992). which has led to the sugges- 
tion that LoBALs are quasars in the process of removing their 
natal cocoons of gas and dust. 

Figure 2 shows a striking trend for more luminous galax- 
ies to have higher outflow velociti es. Similar trends have 
been noted pr eviously. In starbursts. IRupke et alJd2005bl) and 
Martin] d2005l) found strong correlations between outflow ve- 
locity and galaxy mass and star formation rate. Our post- 
starburst galaxies have extraordinarily high outflow velocities 
when compared to their natural analogs, starburst-powered 
ULIRGs and LBGs. Their outflow velocities are compara- 
ble to some of the AGN composite ULIRGs, and at the lower 
end of the range observed for LoBAL quasars. This result im- 
plies that our post-starburst galaxies may harbor both fading 
starbursts and fading quasars. Evidence for the presence of 
an AGN can also be found in the spectra. In massive metal- 
rich galaxies the narrow emission line [O III] A5007 is rela- 
tively uncontaminated by star for mation and a good trac er of 
the AGN's bolometric luminosity dHeckman e t al. 2005b. We 
are able to detect [O III] lines in four of our galaxies using 
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FIG. 2. — Outflow velocity versus absolute B-band magnitude. Our post- 
starburst galaxies (black plus signs) have velocities intermediate between lu- 
minous starbursts and LoBAL quasars. References for the various samples 
are given in the text. 

the SDSS spectra. The galaxies have EW[oni] = 6-9 A and 
L[oni] = 0.4-4 x 10 8 Lq, placing them in the regime of power- 
ful AGN. One of the four also shows [Ne V] A3426 emission, 
which is an unequivocal signpost of AGN activity. 

We achieve better continuum fits to six of our galaxies 
(three of which have [O III] emission) when we add a feature- 
less power-law component {Fx oc A QA ) w ith a spectral slope 
of a\ = -1.6, which is typical of quasars (Vanden Ber k et all 
2001). Without the power-law, high present-day star forma- 
tion rates are implied, which is at odds with the lack of strong 
Balmer emission. The very blue UV continua of these galax- 
ies rules out high dust attenuation as a means of quenching the 
nebular lines. In our best-fit models the power-law supplies 
40-70% of the flux at 3000 A. Curiously, despite the moder- 
ate dilution of the AGN continuum by the galaxy, few spectral 
features characteristic of Type 1 quasars are detected. We are 
able to rule out strong broad H/3 emission lines on the basis 
of the SDSS spectra. Broad Mg II (FWHM-8000 km s" 1 ) is 
present at a low level in SDSS J214000.49+120914.6, but ab- 
sent in the other spectra. The physical reason for the lack of 
broad lines in our AGN-post-starburst composites is unclear. 

Mechanical energy from AGN radio jets has been sug- 
gested as a power source for large-scale outflows (e.g., 
Nesvadba et al. 2006|). For 13/14 galaxies, radio data are 
available from the F aint Images of the Radio Sky at Twenty 
cm survey (FIRST: [Becker etaLl! 1995b . Two galaxies are de- 
tected (Fi 4q Hz = 12, 20 mJy) down to a limit of ~ 1 mJy. 
These two sources show no optical signs of AGN activity, but 
have sufficient radio pow er to be classed as radio-loud AGN 
dKellermann et alj 1 1989t) . Neither galaxy hosts an outflow, 
thus, preliminary evidence disfavors radio jets as the driving 
mechanism. 

It is interesting to consider whether the outflows we observe 
could have been the direct cause of the abrupt drop in the star 
formation rate of our galaxies a few 100 Myr ago. This seems 
plausible if the winds entrained a large fraction of the cold 
ISM. We can crudely estimate the mass in the outflow using 
the Mg II column density that we derive. Because the error 
bars on some individual measurements are large, we use the 
median value, N(Mg + ) = 8.1 x 10 14 cm -2 , in our calculations. 
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The presence of a weak Mg I A2852 line in some of the spec- 
tra implies an ionization correction of ~3%. We translate our 
Mg column into a total gas co lumn by accounting for d eple- 
tion onto dust grains (X=-1.4; Sav age & Sembachll 19961) and 
assuming a Mg/H ratio 2.5 times the solar value. In this way 
we infer N(H)=2 x 10 20 cm -2 , which i s consistent with th e 
median value found for high-z ULIRGs dRupke et al.ll2005bl) . 

The total mass in the outflow depends strongly on how 
far away the absorb ing gas is from the galaxy. Following 
RupkeetalJ d2005blh we assume that the wind is a shell-like 
structure that covers 80% of the optically luminous galaxy. 
To calculate the outer radius of the shell, r out , we adopt a sim- 
ple picture where the wind is launched at the peak of the star 
formation and AGN activity and moves at constant velocity. 
Using the median values of the burst age (f ~ 100 Myr) and 
outflow velocity (v~ 1000 km s" 1 ), we find r oul = vt= 100 kpc. 
The shell's thickness depends on the duration of the outflow. 
We consider two limiting cases: a thick shell with an inner 
radius r m = 1 kpc, and a thin shell with r,„ = 99 kpc. Rupke's 
equation (4) yields masses of M M ,,„^ = 10 9 and 10 11 M© for 
the thick and thin cases respectively, implying that between 1 
and 50% of the galaxies' baryons are in the outflow. The wind 
mass estimated for the thin shell is probably too large: simu- 
lations of starburst and AGN feedback require highly efficient 
and energetic winds to unbind more than 25% of the galaxy's 
initial gas mass (Cox et al., in prep). 

Another possibility we must consider is that the Mg II ab- 
sorber is local to the AGN and does not extend to kilopar- 
sec scales. In this case the wind is unlikely to have played 
a significant role in regulating star formation. For six of our 
galaxies this possibility cannot be ruled out because our es- 
timate of the Mg II covering factor is less than or equal to 
the amount of continuum light contributed by the AGN. We 
obtained data with 5x higher spectral resolution for SDSS 
J082638.41+430529.5in order to obtain a better measurement 
of the covering factor. In the high resolution spectrum, the 
Mg II absorption is near-black at line center, implying that the 
absorber covers both the AGN and the stars. We also detect 
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strong Mg II absorption in three galaxies with no contribution 
to the continuum from an AGN. Hence, in 4/10 galaxies we 
can confirm that the winds are galaxy-scale features indica- 
tive of energetically significant feedback events. Constraints 
on the remaining galaxies await higher resolution spectra. 

6. CONCLUSIONS 

To test currently popular models of early-type galaxy evo- 
lution that incorporate feedback from AGN, we have looked 
for the presence of galactic winds in a sample of massive 
post-starburst galaxies at z = 0.5-0.75. We detect interstel- 
lar Mg II which is blueshifted by 500 - 2000 km s"'in 10/14 
galaxies. These outflow velocities are intermediate between 
those of luminous starbursts and LoBAL quasars, which sug- 
gests that feedback from an AGN may have played a role in 
powering the outflow. In 4/10 galaxies we can confirm that 
the outflows are energetically significant galaxy-wide events, 
and not phenomena local to the AGN. We estimate that the 
outflows reach distances of ^100 kpc and contain upwards 
of 10 9 M Q of gas. We conclude that AGN are likely to have 
played a major role in causing the abrupt truncation of star 
formation in these massive galaxies. 
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TABLE 1 



SDSS Galaxy Name 


z 


M B 


logM, 


Mg II EW 




^max 






(mag) 


(M Q ) 


(A) 


(kms- 1 ) 


(kms" 1 ) 


J081 150.09+471615.2 


0.515 


-22.5 


11.0 


2.5 


-918 ± 27 


-918 ± 27 


J082638.41+430529.5 


0.603 


-23.3 


10.9 


4.8 


-1040 ± 42 


-1232 ± 09 


J082733. 88+29545 1.3 


0.681 


-23.1 


11.3 








J094417. 85+093019.4 


0.514 


-22.7 


10.8 


7.9 


-1245 ± 10 


-1807 ± 13 


J103906.97+453754.1 


0.634 


-23.3 


11.1 


4.1 


-904 ±31 


-1342 ±40 


Jl 10437.46+594639.6 


0.573 


-22.9 


10.9 


3.4 


-1197 ±35 


-1335 ± 22 


Jl 12518.90-014532.5 


0.519 


-23.3 


11.1 


1.9 


-2022 ± 10 


-2022 ± 10 


Jl 14257.23+60371 1.2 


0.568 


-23.6 


11.5 








J124807.16+060111.8 


0.632 


-23.2 


11.2 


2.9 


-489 ± 18 


-489 ± 18 


J150636.30+540220.9 


0.608 


-23.5 


10.9 


2.5 


-1114 ± 66 


-1135 ±78 


J160413.25+393931.4 


0.564 


-23.6 


11.7 








J163541. 72+470924.5 


0.699 


-23.7 


11.5 








J171300.39+281708.2 


0.577 


-23.0 


11.2 


0.8 


-828 ± 35 


-828 ± 35 


J214000.49+120914.6 


0.751 


-23.6 


11.2 


10.4 


-234 ± 40 


-573 ± 48 



NOTE. — We assume Q M = 0.3, = 0.7, and Ho = 70 km s 1 Mpc 1 . M B is ^-corrected to z = and on the Vega system. The Mg II EWs are for the 
interstellar component and are measured in the rest frame. The velocities \ max and v aV g are defined in §4. 



